pdfpiw.uspto.gov/.piw?docid=063771 1 6 

(12) United States Patent 

Mattsson et al. 



US006377116B1 


(io) Patent No.: US 6,377,116 B1 

( 45 ) Date of Patent: Apr. 23, 2002 


(54) PRE-DISTORTER AND CORRESPONDING 
METHOD FOR DERIVING SAME 

(75) Inventors: S. Anders Mattsson, Maineville, OH 
(US); Steve F. Russell, Ames, IA (US) 

(73) Assignee: Iowa State University Research 

Foundation, Inc., Ames, IA (US) 

( * ) Notice: Subject to any disclaimer, the term of this 

patent is extended or adjusted under 35 
U.S.C. 154(b) by 0 days. 

(21) Appl. No.: 09/566,642 

(22) Filed: May 8, 2000 

(51) Int. Cl. 7 H03F 1/26 

(52) U.S. Cl 330/2; 330/149 

(58) Field of Search 330/2, 149; 375/296, 

375/297 

(56) References Cited 

U.S. PATENT DOCUMENTS 

5,867,065 A * 2/1999 Leyendecker 330/149 

5,923,712 A * 7/1999 Leyendecker et al 375/297 

6,118,335 A * 9/2000 Nielson et al 330/2 

6,141,390 A * 10/2000 Cova 375/297 

OTHER PUBLICATIONS 

Maria-Gabriella Di Benedetto and Paolo Mandarini, "Anew 
analog predistortion criterion with application to high effi- 
ciency digital radio links,” IEEE Transactions on Commu- 
nications, vol.43, No. 12, pp. 2966-2974, Dec. 1995. 
Giovanni Lazzarin, Silvano Pupolin and Augusto Sarti, 
"Nonlinearity compensation in digital radio systems,” IEEE 
Transactions on Communications, vol. 42, No. 2/3/4, pp. 
988-999, Feb./Mar./Apr. 1994. 


Aldo N. D ’Andrea, Vincenzo Lottici and Ruggero Reggian- 
nini, “RF power amplifier linearization through amplitude 
and phase predistortion,” IEEE Transactions on Communi- 
cations, vol. 44, No. 11, pp. 1477-1484, Nov. 1996. 

(List continued on next page.) 

Primary Examiner — Steven J. Mottola 

(74) Attorney, Agent, or Firm — Flehr Hohbach Test 

Albritton & Herbert; R. Michael Ananian, Esq. 

(57) ABSTRACT 

A method of deriving a pre-distorter for a corresponding 
non-linear amplifier is disclosed. The first step of the method 
is to form the pre-distorter’s topology. This is done by first 
determining a corresponding inverse topology to the non- 
linear amplifier’s topology. The non-linear amplifier has a 
topology comprising filter and non-linear amplifier ele- 
ments. The inverse topology comprises corresponding 
inverse filter and amplifier elements to the amplifier’s filter 
and non-linear amplifier elements. Then, those of the inverse 
filter and amplifier elements for which compensation with 
the pre -distorter is desired are selected. Any of the inverse 
filter and amplifier elements that were not selected are then 
removed from the inverse topology to form the pre- 
distorter’s topology. Thus, the pre-distorter’s topology com- 
prises the selected inverse filter and amplifier elements. 
After the pre-distorter’s topology is determined, a pre- 
defined input signal is applied to the non-linear amplifier. In 
response, measurements of selected signals of the non-linear 
amplifier are made. Once the measurements are made, 
coefficients that define the selected inverse filter and ampli- 
fier elements’ are estimated. This is done in response to the 
measurements that were made. 
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PRE-DISTORTER AND CORRESPONDING 
METHOD FOR DERIVING SAME 

TECHNICAL FIELD OF THE INVENTION 

The present invention relates generally to pre-distortion of 
an input signal to a non-linear amplifier. In particular, it 
pertains to a pre-distorter for a non-linear amplifier that does 
not actually implement the kernels of a Volterra operator 
used to describe the non-linear amplifier nor the inverse 
kernel of the first kernel of the Volterra operator. 

BACKGROUND OF THE INVENTION 

When designing amplifier, the desired goal is to obtain an 
amplifier without any distortion. This is practically impos- 
sible to do because some of the components, such as 
transistors, tubes, etc., used in an amplifier are non-linear. 
Furthermore, the capacitors and inductors used in an ampli- 
fier cause frequency dependent distortion. For an RF (i.e., 
radio frequency) amplifier, this results in two different 
effects. One effect is that the output signal is distorted at the 
carrier frequency and the other effect is that there are 
harmonics in the output signal at multiples of the carrier 
frequency. In order to deal with such distortion, the input 
signal can be pre-distorted (i.e., pre-shaped) with a pre- 
distorter before it is sent to the nonlinear amplifier. 

To derive a corresponding pre-distorter for a non-linear 
amplifier, a standard approach is to describe the amplifier in 
terms of a Volterra operator . The advantage to doing so is 
that, given that the inverse kernel of the first kernel of the 
Volterra operator for the non-linear amplifier exists, the 
corresponding pre-distorter can be derived by describing it 
in terms of yet another Volterra operator. A good discussion 
of Volterra operators can be found in Schetzen, Martin, The 
Volterra and Wiener Theories Of Nonlinear Systems, John 
Wiley & Sons, New York, 1980, which is incorporated by 
reference herein. 

The Volterra operator for a pre-distorter derived using this 
approach will be in terms of the kernels of the Volterra 
operator for the corresponding non-linear amplifier and the 
inverse kernel of the first kernel of the Volterra operator for 
the nonlinear amplifier. Unfortunately, this means that a 
pre-distorter will have to implement these kernels. As evi- 
denced in The Volterra and Wiener Theories Of Nonlinear 
Systems, this is very messy and impractical since determin- 
ing the kernels is in general awkward and, as the order 
increases, the practical implementation of the kernels 
becomes excessively complicated. Therefore, there is a need 
for an efficient method for deriving pre-distorters, particu- 
larly those with high orders, that do not actually implement 
these kernels. 

SUMMARY OF THE INVENTION 

In summary, the present invention comprises a method of 
deriving a pre-distorter for a corresponding non-linear 
amplifier. The present invention also comprises the pre- 
distorter itself. 

The first step of the method is to define the pre-distorter’s 
topology. This is done by first determining a corresponding 
inverse topology to the non-linear amplifier’s topology. The 
non-linear amplifier has a topology comprising filter and 
non-linear amplifier elements. The inverse topology com- 
prises corresponding inverse filter and amplifier elements to 
the amplifier’s filter and non-linear amplifier elements. 
Then, each of the filter and non-linear amplifier elements of 
the amplifier for which compensation with the pre-distorter 
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is desired are identified. Any of the inverse filter and 
amplifier elements corresponding to one of the filter and 
non-linear amplifier elements that is not identified is then 
removed from the inverse topology. This forms the pre- 
5 distorter’s topology so that it comprises the corresponding 
inverse filter and amplifier elements to the identified filter 
and non-linear elements. 

The next step of the method is to define the selected 
inverse filter and amplifier elements’ responses. To do so, an 
10 input signal is applied to the non-linear amplifier. In 
response, measurements of selected signals of the nonlinear 
amplifier are made. Once the measurements are made, 
coefficients that define the selected inverse filter and ampli- 
fier elements’ responses are estimated. This is done in 
15 response to the measurements that were made. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIGS. 1 to 3 show block diagrams of amplifier systems 
20 comprising non-linear amplifiers and corresponding pre- 
distorters in accordance with the present invention. 

FIG. 4 shows a block diagram of the configuration of the 
amplifier of FIG. 1. 

FIG. 5 shows a block diagram of the configuration of the 
25 pre -distorter of FIG. 1. 

FIG. 6 shows the performance of the pre-distorter of FIG. 

5. 

FIG. 7 shows a block diagram of the configuration of the 
30 amplifier of FIG. 2. 

FIG. 8 shows a block diagram of the configuration of the 
pre-distorter of FIG. 2. 

FIG. 9 shows a block diagram of the configuration of the 

amplifier of FIG. 3. 

35 

FIG. 10 shows a block diagram of the inverse configu- 
ration of the amplifier of FIG. 3. 

FIG. 11 shows a block diagram of the configuration of the 
pre-distorter of FIG. 3. 

40 FIG. 12 shows the performance of the pre-distorter of 
FIG. 11. 

DETAILED DESCRIPTION OF THE 
INVENTION 

45 

The invention disclosed herein provides a novel method 
for deriving pre-distorters for a broad class of non-linear 
amplifiers. To illustrate the method, three specific examples 
are provided. Referring to FIGS. 1 to 3, in these examples, 
50 pre-distorters 100, 200, and 300 are derived for respectively 
a simple non-linear amplifier 101, a multi-cascaded non- 
linear amplifier 201, and an EER (i.e., envelope elimination 
and restoration) non-linear amplifier 301. 

55 Pre-Distorter 100 for Simple Non-linear Amplifier 

101 

FIG. 4 shows the topology of the simple non-linear 
amplifier 101. The amplifier 101 comprises a non-linear 
60 amplifier element 102 without memory and a linear filter 
element 103 with memory. The amplifier element 102 and 
the linear filter element 103 are cascaded (i.e., connected in 
series). 

The amplifier element 102 has a characteristic non-linear 
65 response f(x(t)) to an input signal x(t). Thus, it is clear that 
the output signal z(t) of the amplifier element 102 in 
response to an input signal x,„(t) is given by: 
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z(t)-f(x,„(t)) (1) 

where x,.„(t)=x(t). 

The filter element 103 has a characteristic linear response 
h(t). This means that the output signal y 4 jv (t) of the filter 
element 103 (and also the amplifier 101) can be expressed 
as: 

Wt)-Kt)*f(*Jt)) (2) 

The output signal y dis (i) may be distorted by the non- jo 
linearity of the amplifier element 102 and/or the bandlimit- 
ing effect of the filter element 103. 

Derivation of Pre-Distorter 100 

The first step in deriving the corresponding pre-distorter 15 
100 in FIG. 5 to the amplifier 101 in FIG. 4 is to define the 
topology of the pre-distorter 100. This is done by first 
determining the inverse topology to that of the amplifier 101. 
This inverse topology comprises a corresponding inverse 
amplifier element 104 to the amplifier element 102 and a 20 
corresponding inverse filter element 105 to the filter element 
103. 

Then, each of the amplifier and filter elements 102 and 
103 of the amplifier 101 of FIG. 4 for which compensation 
is desired are identified. Each inverse element 104 or 105 25 
corresponding to an element 102 or 103 for which compen- 
sation is not desired is then removed to form the topology of 
the pre-distorter 100. 

In the example provided here, the amplifier element 104 ^ 
is assumed to have a significant non-linear behavior at the 
carrier frequency of the input and output signals x,„(t) and 
y 4 is (t) and the filter element 103 is assumed to be a bandpass 
filter with a significant bandlimiting effect at the carrier 
frequency. In this case, the distorted output signal y dis (t) is 
significantly distorted at the carrier frequency by the non- 
linear behavior of the amplifier element 102 and the band- 
limiting effect of the filter element 103. As a result, com- 
pensation is desired for both the amplifier and filter elements 
102 and 103. As shown in FIG. 5, both the inverse amplifier 
and filter elements 104 and 105 are therefore used in the 
topology for the pre-distorter 100. 

Still referring to FIG. 5, to compensate for the response 
h(t) of the filter element 103, the corresponding inverse filter 
element 105 has a characteristic response h - 1 (t). The output 45 
signal w(t) of the filter element 105 in response to an input 
signal x 1 >f (t) is therefore given by: 

w(t)=h- 1 (t)*x,.„(t) (3) 

The inverse amplifier element 104 has a characteristic so 
response f _ 1 (x(t)) to an input signal (x(t)). This is used to 
compensate for the response f(x(t)) of the corresponding 
amplifier element 102. In view of Eq. (3), the pre-distorted 
output signal x pK Jf) of the inverse filter element 105 (and 
also the pre-distorter 100) in response to the output signal 55 
w(t) of the filter element 105 is given by: 


N 


~ Yj b " x »> W" 

n = 0 


(5) 


where the set of coefficients {b 0 , . . . .b^} defines the 
response f(x(t)) of the amplifier element 102. Assuming that 
the inverse response f _ 1 (f(x,-„(t))/g) to the response f(x irt (t)) 
exists and is bounded by an input signal f(x ;> 1 (t))/g, the 
behavior of the amplifier element 102 can also be modeled 
in an inverted manner with the following series expansion: 


Xin(t) = f \f(x m (t)) 



f{Xm(t'l) 

s 


( 6 ) 


where P is the selected order of the pre-distorter 100 and the 
set of coefficients {d 0 , . . . ,d P } defines the response D 1 (x(t)) 
of the inverse amplifier element 104 of FIG. 5. The gain g 
is the desired gain of the amplifier element 102 of FIG. 4. It 
appears in Eq. ( 6 ) so that the coefficients {d 0 , . . . ,d p } of the 
response D 1 (x(t)) are properly defined not to cause ampli- 
fication of the pre-distorted output signal x pre Jt) by the gain 
g- 

If the filter element 103 is assumed to be minimum phase, 
its behavior can be modeled by an AR (i.e., auto regressive) 
filter. Thus, the corresponding Z transform H(z) for the 
response h(t) can be modeled by the series expansion: 


M 

£ a m z~ m 

m = 0 


where a 0 =l and the set of coefficients {a 0 , . . . ,a M } defines 
the response h(t). Furthermore, the fact that the filter element 
103 is assumed to be minimum phase means that the inverse 
response h _ 1 (t) to the response h(t) exists. Thus, this 
response h - 1 (t) of the inverse filter element 105 of the 
pre-distorter 100 of FIG. 5 is also defined by the coefficients 
{a 0 , . . . ,a M } since its Z transform H _ 1 (z) can be modeled 
by the series expansion: 


M 

tf(z)- 1 

m~ 0 


( 8 ) 


As just mentioned, the response h J (t) does exists. Thus, 
Eq. (2) can be re-written as follows: 

ya i ,(t)*h- 1 (t)=f(x,„(t)) (9) 

In view of Eqs. (5) and ( 8 ), Eq. (9) can be re-written to 
model the combined behavior of the amplifier and filter 
elements 102 and 103 of the amplifier 101 of FIG. 4 in the 
following manner: 


W^UwttKMoxtt)) (4) 

where w(t)=x(t). 

The next step in deriving the pre-distorter 100 is to define 
the responses f _ 1 (x(t)) and h - 1 (t) for the inverse amplifier 
and filter elements 104 and 105 of FIG. 5. In order to do so, 
it is necessary to develop a model of the behavior of the 
corresponding amplifier and filter elements 102 and 103 of 
FIG. 4. This may be done in the following manner. 

The behavior of the amplifier element 102 can be modeled 
(i.e., approximated) by the series expansion: 


N M (10) 

ydis(t) ~ ^ b n Xj n (t) n - ^ a m yss<,t - rnT) 

n = 0 m = 1 

where 1/T represents a sampling rate. 

Using the assumption that the filter element 103 is mini- 
mum phase, the model of the amplifier 101 represented by 
Eqs. ( 6 ) and (10) is linear in its parameters. Thus, it can be 
easily and efficiently estimated using an LMS (i.e., least 
mean squares) estimation. 
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In performing this LMS estimation, an input matrix Xj 
and an output vector Y must first be defined as follows: 


6 

However, in order to define the response f 1 (x(t)) of the 
inverse filter element 104 of FIG. 5, an additional LMS 


xM° 

xM 1 

- x,„(tf 

y disit -T) ■ 

■ ya.lt -MT) 

(ID 

X; n {t+T)° 

X m u + T)' ■ 

- x in (t+Tf 

ya.lt -T + T) ■ 

• ya.lt -MT + T) 


x in (t + 2Tf 

ydisW 
yssif+T) 
y<Ss(t + 2T) 

*,„(; + 27y ■ 

■■ Xi„(t + 2T) N 

ya.lt -T + 2T) ■ 

• ya.lt -MT + 2T) 

(12) 


As will be discussed later in more detail, the elements of the 
matrix X-j and the vector Y are obtained by applying an input 
signal x in (t) to the amplifier 101 . This input signal x,„(t) 
should be the same type of input signal (i.e., have the same 20 
bandwidth and amplitude distribution) that will be applied to 
the amplifier system 106 of FIG. 1 during actual operating 
conditions. This will ensure that the responses l^ 1 (x(t)) and 
h -1 (t) for the inverse amplifier and filter elements 104 and 
105 of FIG. 5 are specifically defined for this type of input 
signal. Measurements of the input signal x,„(t) and the 
output signal y dj> (t) are then made from which the elements 
of the matrix Xj and vector Y are obtained. These measure- 
ments are made at the sampling rate 1/T mentioned earlier. 

Based on Eq. (10) of the model, the input matrix X x and 
the output vector Y are related to the coefficient vector 30 
[b 0 . . . b^-a-L . . . -a M Y by the following expression: 


estimation must be made. In performing this additional LMS 
estimation, the output signal z(t) of the amplifier element 
104 of the amplifier 101 of FIG. 4 must be estimated in 
response to the input signal x,„(t). This is done by defining 
the corresponding output vector Z t : 


z(r) 

z(r + T) 
Z(t + 2T) 


and then computing it according to: 


bo 


(17) 


y = x, 


bo 

b N 


—ai 


(13! 
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Zi = X, 


b N 

0 


0 


An output matrix Z 2 and an input vector X 2 can then be 
40 defined as follows for this LMS estimation: 


In order to obtain the coefficient vector [b 0 . . . b^aj . . . 
-a m Y the LMS solution for Eq. (13) can be easily computed 
according to: 


(14) 


45 


b N 

-ai 


= (x?Xi) ‘x^r 


~ a M 


50 


z(r)°/g° z(t)'lg l 
z(t+Tf/g 0 z(r + T) l /g l 
z(r + 2Tf/g° zlt + 2T) 1 /g 1 


z(tf/g p 

Z(t + T) p lg p 
z(r + 2Tf/g p 


(18) 


Xin (0 
Xinit+T) 

x in (t + 2T) 


(19) 


where the matrix X^ is the Hermitian transpose of the 
matrix X 2 . 

M (15) 

H(z) x^a m z~ m 

m=0 


The elements for the matrix Z 2 and the vector X 2 are 
obtained from the elements of the vector Z 2 and the matrix 
55 X 1? respectively. 

Based on Eq. (6) of the model, the output matrix Z 2 and 
the input vector X 2 are related to a coefficient vector [d 0 . . . 
d according to the following expression: 


In this manner, the sets of coefficients {b 0 , . . . jb^J and 
{a 0 , . . . ,a M } are estimated for properly defining the 
responses f(x(t)) and h(t) of the amplifier and filter elements 
102 and 103 of the amplifier of FIG. 4. The response h -1 (t) 
for the inverse filter element 105 of the pre-distorter 100 of 65 
FIG. 5 is therefore also defined by the 5 coefficients 
{a 0 , . . . ,a M }. 



( 20 ) 


In order to obtain the coefficient vector [d 0 . . . d P ] r , the LMS 
solution for Eq. (20) can be computed using matrix inversion 
according to: 
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-- (z%z 2 y l z%x 2 


(21) 


5 


where Z 2 is the Hermitian transpose of Z 2 . Therefore, in this 
approach, the set of coefficients {d 0 , . . . ,d P } are directly 
estimated. io 


In this manner, the set of coefficients {d 0 , . . . ,d P } is 
estimated for properly defining the response f" 1 (x(t)) of the 
inverse amplifier element 105 . Since the response h -1 (t) of 
the inverse filter element 104 was previously defined with 
the estimated set of coefficients {a 0 , . . . ,a M }, the behavior I s 
of the inverse amplifier and filter elements 104 and 105 of 
the pre-distorter 100 can now be properly modeled by two 
cascaded processes. The first process is performed by the 
inverse filter element 105 and can be expressed as follows: 


w(r) * [ x-m(t - T ) x in (t - 2 T) 


-a. 




-«2 


(22) 


-am 


25 


The second process is performed by the inverse amplifier 
element 104 and can be expressed in the following manner: 


XpndW* N0° MO* M0 P ] 


( 23 ) 


30 


35 


Since these processes only involve simple additions and 
multiplications, the implementation of the pre-distorter 100 
with the inverse filter and amplifier elements 104 and 105 is 
also simple and efficient. 


Amplifier System 106 and Operation 

FIG. 1 shows an amplifier system 106 that comprises the 
amplifier 101 and the corresponding pre-distorter 100 that 
was just derived. The pre-distorter 100 pre-distorts the 45 
original input signal x ;>1 (t) according to Eqs. (21) and (22) to 
produce the pre-distorted input signal x pre Jt). The pre- 
distorted input signal x pred (t) is then amplified by the ampli- 
fier 101 to produce the undistorted output signal y 0 „ ; (t). The 
output signal y our (t) is not distorted by the non-linear behav- 50 
ior of the amplifier element 102 and the bandlimiting effect 
of the filter element 103 that were compensated for with the 
pre-distorter 100 . 

Performance Evaluation of Pre-Distorter 100 
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-continued 


h(x(t)) = 


2 x{t) 

1 + x{t ) 2 


iicxjl) 

c 3(1+j(i)2) 


(25) 


The first response f 2 (x(t)) models an arbitrary amplifier. The 
second response f 2 (x(t)) models a traveling wave tube. This 
response is complex valued and models the AM to FM (i.e., 
amplitude modulation to frequency modulation) conversion 
that takes place in the amplifier element 102 . 

For both output responses f 2 (x(t)) and f 2 (x(t)) used for the 
amplifier element 102 , the same characteristic response h(t) 
was used for the filter element 103 of the amplifier 101 . The 
Z transform H(z) for this response h(t) is expressed by: 


H(l) = 


1 

1 + (- 1.7798 + 0 . 1176/);- 1 + 
( 1.3885 - 0 . 1046 /) z " 2 + 


(- 0.5340 + 0 . 0353 /);-' + 0 . 090;" 4 


( 26 ) 


The response h(t) is complex valued and is used as the 
baseband equivalent of an RF filter. It causes significant 
distortion in the output signal y Ai (t). 

In the simulations, two input signals x 2 (t) and x 2 (t) were 
used as the input signal x ;>I (t). The first input signal x 2 (t) 
comprised two independent bandlimited white Gaussian 
noise components. One component was used as the I com- 
ponent (i.e., the real component) of the signal x 2 (t) and the 
other component was used as the Q component (i.e., com- 
plex component) of the signal x 2 (t). The second input signal 
x 2 (t) comprised a normalized two-tone single-sideband sig- 
nal. The fact that the response f 2 (x(t)) is not invertible for all 
inputs does not cause any problems in the estimation process 
due to the way in which the estimations are made. However, 
both signals Xj(t) and x 2 (t) were normalized to have unity 
amplitude since the response f 2 (x(t)) is only invertible in the 
range 0 to 1. 

For each of the four combinations of responses f 2 (x) and 
f 2 (x) and input signals x 2 (t) and x 2 (t), a corresponding 
simulation of the amplifier 100 was performed. Each of 
these simulations also used the response h(t) just described. 
For each simulation, measurements of the input signal x 2 (t) 
or x 2 (t) and the corresponding distorted output signal y rfij (t) 
were made for obtaining the elements of the matrices and 
vectors X 2 , Y, Zi, Z 2 , and X 2 . Then, for each simulation, 
pre-distorters 100 with orders P=3, 5, 7, 9, and 11 were 
derived in the manner described earlier. 

To determine the performance of the pre-distorters 100 
derived for each simulation, an SDR (signal to distortion 
ratio) was used. The SDR can be defined as: 


SDR = - lOloglO 


E(ydi,(t)-gxjt)) 2 ' 

E(gxjr)) 2 


E(x(t)f = 



( 11 ) 


( 28 ) 


The method just described for deriving the pre-distorter 
100 was verified with simulations using two different char- 
acteristic responses fj(x(t)) and f 2 (x(t)) for the response 
f(x(t)) of the non-linear amplifier element 102 of the ampli- 
fier 101 of FIG. 4. 

These responses f 2 (x(t)) and f 2 (x(t)) are given by: 


where E(x(t)) 2 is simply the energy of a signal x(t) and is 
given by: 

FIG. 6 shows the performances. For each simulation, the 
improvement in SDR when going from a pre-distorter 100 of 
order P=3 to order P=ll is a minimum of 50 dB. This 
demonstrates that much improved pre-distorters 100 can be 
achieved when increasing the order P. 


60 


.v(r) 

= X ( t ) + - j ^ j - 


( 24 ) 65 Alternative Estimation Techniques 

The set of coefficients {d 0 , . . . ,d P } were directly 
estimated using Eq. (21) in the example provided earlier. 
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However, the set of coefficients {d 0 , . . . ,d P } can be 
indirectly estimated as well using a similar approach for the 
LMS estimation. In this approach, an output matrix Z 3 can 
be defined as follows: 

zU)° zU) l> - z(t) p (29) 

Ztt + T? z(t + T) 1 ■■■ z(t+Tf 

Z 3 = 

z(r + 27') 0 z(r + 2r)‘ ■■■ z(r + 2 T) p 


where the elements of the matrix Z 2 are not adjusted for the 
gain g. As with the elements of the matrix Z 2 , the elements 
of the matrix Z 3 are obtained from the elements of the vector 
Z 3 . 

The output matrix Z 3 and the input vector X 2 are related 
to a coefficient vector [c 0 . . . Cp] 7 as follows: 


co 


X 2 =Z 3 


cp 


(30) 


The LMS solution for Eq. (30) is then computed according 
to: 


c 0 (31) 

: =(z('z 3 r 1 zj'x 2 

Cp. 


where Z 3 H is the Hermitian transpose of Z 3 . After the 
coefficient vector [c 0 . . . c P ] 7 is obtained from the LMS 
solution, the coefficient vector [d 0 . . . dp] 7 is obtained 
according to: 


do 


' Colg 0 ' 

di 

= 

Cl/g l 

d P 


. Cp/g p 


where dj=Cj/g J . In this way, the set of coefficients {d 0 , . . . 
,d P } are indirectly estimated. 

Furthermore, in the LMS estimations discussed earlier for 
obtaining the coefficient vectors [b 0 . . . b w -a 3 . . . -a M ] , 
[d 0 . . . d p ] 7 , and [c 0 . . . Cp] 7 , matrix inversions were used 
to solve Eqs. (14), (20), and (30). But, these matrix inver- 
sions can be avoided by iteratively finding these coefficient 
vectors [b 0 . . . b^-a, . . . -a M ] 7 , [d 0 . . . dp] 7 , and [c 0 . . 
. Cp] 7 , using RLS (i.e., recursive least squares) estimations 
or GAL (i.e., gradient adaptive lattice) estimations in the 
manner described in Haykin, Simon , Adaptive Filter Theory, 
2nd. Ed. Prentice Hall, Englewood Cliffs, N.J., 1991, which 
is hereby incorporated by reference. Furthermore, SVDs 
(i.e., singular value decomposition) of the matrices X-^Xj, 
Z^Z-,, and Z 3 77 Z 3 in Eqs. (14), (20), and (30) can be used. 
This becomes useful when the matrices X/LX.,, Z^Z z , and 
Z 3 h Z are close to singular. In practice, this will happen 
sooner or later as the selected order P of the model of the 
amplifier 101 (and the derived pre-distorter) increases. But, 
generally this is an indication that the order P is too high. 

Alternative Configurations for Pre-distorter 100 

In the earlier example, compensation was desired for both 
the amplifier and filter elements 103 and 104. Therefore, the 
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corresponding inverse amplifier and filer elements 104 and 
105 were utilized in the topology of the pre-distorter 100. 

However, as those skilled in the art will recognize, it may 
be desired to compensate for only the amplifier element 103 
5 or only the filter element 104. For example, if only the 
non-linear behavior of the amplifier element 103 or only the 
bandlimiting effect of the filter element 104 causes signifi- 
cant distortion in the output signal y Ai (t) at the carrier 
frequency, then only the corresponding inverse amplifier 
10 element 104 or inverse filer element 105 would be used in 
the topology of the pre-distorter 100. The corresponding 
responses f _1 ((x(t)) or h -1 (t) would still be derived in the 
same manner as described previously. But, only the distor- 
tion of the amplifier or filter element 102 or 103 that was 
I 5 compensated for with the pre-distorter 100 would not distort 
the output signal y ouf (t). 

Pre-distorter for Multi-Cascaded Non-Linear 
Amplifier 

20 Referring to FIG. 7, the method just described for deriv- 
ing the pre -distorter 100 of FIGS. 1 and 5 can be applied to 
a multi-cascaded non-linear amplifier 201 with cascaded 
non-linear sub-amplifiers 101 3 to 101 7 , where the index i=l 
to I. Each sub-amplifier 101,- is configured and operates like 
25 the amplifier 101 of FIGS. 1 and 4 because it comprises a 
corresponding non-linear amplifier element 102,- without 
memory and a corresponding linear filter element 103, with 
memory. 

30 Derivation of Pre-Distorter 200 

Like the derivation of the pre-distorter 100, the first step 
in deriving the pre-distorter 200 in FIG. 8 is to define its 
topology. The topology of the pre-distorter 200 is formed by 
first determining the inverse topology to that of the amplifier 
201. As a result, the inverse topology comprises cascaded 
sub-pre-distorters 100 7 to 100,, where there is a correspond- 
ing sub-pre-distorter 100, for each sub-amplifier 101,-. Each 
sub-pre-distorter 100 ,- is configured, derived, and operates 
4Q like the pre-distorter 100 of FIGS. 1 and 5 because it 
comprises a corresponding inverse amplifier element 104,- to 
the amplifier element 102, and a corresponding inverse filter 
element 105, to the filter element 103,-. 

To further define the topology of the pre-distorter 200, 
45 each element 103/ and 104/ of a sub-amplifier 101/ for which 
compensation is desired is identified. Then, any inverse 
element 104, or 105,- corresponding to an element 103,- or 

104, of a sub-amplifier 101, for which compensation is not 
desired is removed from the inverse topology. This includes 

50 removing a sub-amplifier 101,- entirely if compensation is 
not desired for either of its inverse amplifier and filter 
elements 104, and 105,. 

In the example provided here, compensation is desired for 
both the amplifier and filter elements 102,- and 103,- of each 
55 sub-amplifier 101,- for the same reasons as those expressed 
earlier for the amplifier and filter elements 102 and 103 of 
the amplifier 101 of FIG. 4. Thus, as shown in FIG. 8, both 
the inverse amplifier and filter elements 104,- and 105,- are 
used in the topology for the corresponding sub-pre-distorter 
60 100,. 

The next step is to define the responses f 31 1 ,(x(t)) and 
h _1 ,-(t) for the inverse amplifier and filter elements 104,- and 

105, of each sub-distorter 100, of FIG. 8. Referring to FIG. 
7, this is done by applying an input signal x,„(t) to the first 

65 sub-amplifier 101 3 of the amplifier 201. Here as well, the 
input signal x,„(t) should be the same type of input signal 
that will be applied to the amplifier system 206 of FIG. 2 
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during actual operating conditions. This will ensure that the 
responses f _1 ,(x(t)) and h _1 ,-(t) for the inverse amplifier and 
filter elements 104, and 105, of each sub-distorter 100, are 
specifically defined for this type of input signal. 

For each sub-amplifier 101,- except for the first sub- 
amplifier 101], the corresponding input signal is the output 
signal y<ft,,_i(t) of the previous sub-amplifier 101,- For 
each sub-amplifier 101,-, measurements of the corresponding 
input signal x,„(t) or y dis ^(t) and the corresponding output 
signal y dis ,.(t) of the sub-amplifier 101,- are then made. Based 
on these measurements, the responses l^ 1 ,-(x(t)) and h -1 ,^) 
for the inverse amplifier and filter elements 104, and 105,- of 
the sub-amplifier 101,- are defined. This is done in the same 
manner in which the responses f _1 (x(t)) and h -1 (t) for the 
inverse amplifier and filter elements 104 and 105 of the 
pre-distorter 100 of FIG. 5 were defined except for a few 
minor differences discussed next. 

For each sub-amplifier 101,, the terms z(t) and y dis (t) are 
replaced by the terms z,(t) and y dis ,(t) to properly identify 
the corresponding output signal z,-(t) and the corresponding 
distorted output signal y djj! (t) of the amplifier and filter 
elements 102,- and 103,, respectively, of the sub-amplifier 
101,, Similarly, the terms f(x(t)) and h(t) are replaced by the 
terms f,-(x,-(t)) and h,-(t) to properly identify the characteristic 
responses f,(x,(t)) and h,-(t) of the amplifier and filter ele- 
ments 102,- and 103,, In addition, for each sub-amplifier 101, 
except for the first sub-amplifier 101,, the term x,„(t) is 
replaced by the term y rflJ ,_ 1 (t) to properly identify the 
corresponding input signal y rf£Sji _i(t) to the amplifier element 
102 ,- of the sub-amplifier 101,, 

Furthermore, for each pre-distorter 100,, the terms w(t) 
and x pre J$) are replaced by the terms w,(t) and x predt {i) to 
properly identify the corresponding internal signal w,(t) and 
the corresponding output signal x pred ,{\) of the inverse filter 
and amplifier elements 105, and 104,, respectively, of the 
sub-pre-distorter 100,, Also, the terms f(x _1 (t)) and h -1 (t) are 
replaced by the terms f(x, _1 (t)) and h, _1 (t) to properly 
identify the characteristic responses f(x, _1 (t)) and h,- -1 (t) of 
the inverse amplifier and filter elements 104, and 105,, 
Finally, for each sub-pre-distorter 100, except for the first 
sub-pre-distorter 100„ the term x,„(t) is replaced by the term 
x pred i+i(t) t0 properly identify the corresponding input sig- 
nal x prerf ,_ 1 (t) to the inverse filter element 105, of the 
sub-pre-distorter 100 ,, 

The sub-pre-distorters 100, to 100] are then cascaded in 
opposite order to that of the sub-amplifiers 101] to 100, to 
form the pre-distorter 200. For each sub-pre-distorter 100, 
except the final sub-pre-distorter 100„ the output signal 
x predA t) i s the input signal x ;> ,,-_ 1 (t) to the next sub-pre- 
distorter 100, The final output signal x ^(t) from the 
final sub-pre-distorter 100, is the pre-distorted input signal 
to be provided to the amplifier 201. 

Amplifier System 206 Topology and Operation 

FIG. 2 shows an amplifier system 206 that comprises the 
amplifier 201 and the corresponding pre-distorter 200. An 
original input signal x,„(t) is applied to the pre-distorter 200 
and pre-distorted by the pre-distorter 200 to produce the 
pre-distorted input signal x pKd ,(t) provided to the amplifier 
201. In doing so, each sub-pre-distorter 100, pre-distorts the 
corresponding input signal x,-„(t) or x pn , di+1 (t) to produce the 
corresponding output signal x predl (\) according to Eqs. (21) 
and (22), but where x p „ du] (t)=xjl)) for i<I, w,-(t)=w(t), 
and x predl {t)=x pri ,Jt). The final output signal x pred l (t) from 
the final sub-pre-distorter 100, is the pre-distorted input 
signal provided to the amplifier 201. 
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The pre-distorted input signal x rfl (t) is then amplified 
by the amplifier 101 to produce the undistorted output signal 
y ««</(!)• In this process, the corresponding input signal for 
the first sub-amplifier 101, is the pre-distorted input signal 
5 x rfl (t) while the corresponding input signal for every 
other sub-amplifier 101,- is the distorted output signal y out 
i-i(t) of the previous sub-amplifier 101, Each sub- 
amplifier 101, amplifies the corresponding input signal x d 
i(t) or y ou ,,-_ 1 (t) to produce the corresponding output signal 
to y oul X t) The final output signal y OM ,j(t) from the final sub- 
amplifier 100, is the undistorted output signal y outd {l) This 
output signal y OI „_Xt) is not distorted by the non-linear 
behavior of any amplifier element 102,- and the bandlimiting 
effect of any filter element 103,- that were compensated for 
15 with the pre-distorter 200. 

Performance Evaluation and Alternative 
Configurations and Estimation Techniques 

20 As mentioned earlier, each sub-pre-distorter 100,- is 
configured, derived, and operates like the pre-distorter 100 
of FIGS. 1 and 5. Thus, the performance evaluation and 
alternative configurations and estimations discussed earlier 
for the pre-distorter 100 apply to each sub-pre-distorter 100,, 

25 

Pre-distorter for EER Non-Linear Amplifier 

FIG. 9 shows the topology of the EER non-linear ampli- 
fier 301. The amplifier 301 comprises an AM removal 
30 element (e.g., a hard limiter) 302, a non-linear amplifier 
element 303 without memory, an envelope detection ele- 
ment 304, a first filter element 305 with memory, and a 
second filter element 306 with memory. The AM removal 
element 302, the amplifier element 303, and the second filter 
35 element 306 are cascaded to form a main path of the 
amplifier 301. The envelope detection element 304 and the 
first filter element 305 are cascaded to form a control path of 
the amplifier 301. In the control path, the first filter element 
305 is connected to the amplifier element 303 for controlling 
40 the amplification by the amplifier element 303 in the main 
path. 

The AM removal element 302 and the envelope detection 
element 304 both receive an input signal x,„(t) and split it 
45 into a PM (pulse modulation) signal x PA ^t) and an envelope 
signal x em ,(t). The envelope signal x e „ v (t) is bandlimited by 
the first filter element 305 to produce a control signal x co „(t) 
for controlling the amplification by the amplifier element 
303. The filter element 305 has a characteristic response 
so h,(t) such that the control signal x con (t) is given by: 

*e»»(t)-6l(0**o.v(t) (33) 

The PM signal x PM (t) is amplified by the amplifier 
55 element 303 in response to the control signal x co „(t) to 
produce the amplified signal y amp (t). The amplifier element 
303 has a characteristic complete response x a (t)f(x fc (t)) to a 
first input signal x a (t) and a second input signal x fc (t) that 
controls amplification of the first input signal x a (t). The 
60 amplified signal y „„,,,( t) in response to the PM signal x PM ( t) 
and the control signal x co „(t) is given by: 

y« mp (b=Xp M (t)f(x c „„(t)) (34) 

65 where W0=x a (t) and x co Ji)=x b ({). 

The complete response x a (t)f(x fc (t)) of the amplifier ele- 
ment 303 includes the characteristic non-linear partial 
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response f(x fc (t)). In view of the partial response f(x fc (t)), the 
envelope signal y cnv (y) for the amplified signal y om „(t) is 
given by: 

y„v(t)=f(x COT (t)) (35) 

The amplified signal y amp (i) and the envelope signal y env (y) 
may be distorted by the non-linear behavior of the amplifier 
element 303 and/or the bandlimiting effect of the filter 
element 305. The amplified signal y arnp (l) is then bandlim- 
ited by the second filter element 306 to produce an output 
signal y As (t) of the amplifier 301. The filter element 306 has 
a characteristic response h 2 (t). Thus, the output signal y rfjv (t) 
can be expressed by: 

Ya,i(0=li2(f>*y» m/ ,(d (36) 

Because the output signal y rf£s (t) is based on the amplified 
signal y amp (t), the output signal y rffc (t) may be distorted by 
the non-linear behavior of the amplifier element 303 and/or 
the bandlimiting effect of the filter element 305. 
Furthermore, the output signal y di Jt) ma y also be distorted 
by the bandlimiting effect of the filter element 306. 

Derivation of Pre-Distorter 300 

As in the derivation of the pre-distorters 100 and 200, the 
first step in deriving the pre-distorter 300 in FIG. 10 is to 
define its topology. This is done by determining the inverse 
topology to that of the corresponding amplifier 301 in FIG. 
9. This inverse topology is shown in FIG. 11 and comprises 
an AM removal element 302 and an envelope detection 
element 304 like those in the amplifier 301. This inverse 
topology also comprises a mixer 307, a corresponding 
inverse amplifier element 308 to the amplifier element 303, 
a corresponding first inverse filter element 309 to the first 
filter element 305, and a corresponding second inverse filter 
element 310 to the second filter element 306. The envelope 
detection element 304, the inverse amplifier element 308, 
and the first inverse filter element 309 are cascaded to form 
a first path in the inverse topology. The AM removal element 
302 forms a second path of the inverse topology. The first 
and second paths are connected to the second inverse filter 
via the envelope detection element 308 and the AM removal 
element 302, respectively. The first and second paths are 
combined with the mixer 307. 

Similar to the derivation of the pre-distorters 100 and 200, 
element 303, 305, and 306 of the amplifier 301 of FIG. 9 for 
which compensation is desired is identified. Then, any 
inverse element 308, 309, or 310 corresponding to an 
element 303, 305, or 306 of the amplifier 301 of FIG. 9 for 
which compensation is not desired is removed from the 
topology to form the topology of the pre-distorter 100. 

In the example presented here, the amplifier element 303 
is assumed to have a significant non-linear behavior at the 
carrier frequency of the input and output signals x iM (t) and 
y di JS). This is due to the partial response f(x t (t)) of the 
amplifier element 301. Furthermore, the first filter element 
305 is assumed to be a bandpass filter whose response h,(t) 
causes a significant bandlimiting effect in the output signal 
y di JS) at the carrier frequency. However, the second filter 
element 306 is assumed to be a bandpass filter whose 
response h 2 (t) suppresses harmonic frequencies in the output 
signal y dis (t) with no significant bandlimiting effect at the 
carrier frequency. This means that only the partial response 
f(x t (t)) and the response h 1 (t) of the amplifier and filter 
elements 303 and 305, respectively, cause significant dis- 
tortion in the distorted output signal y dlv (t) at the carrier 
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frequency. Compensation is therefore desired for the ampli- 
fier and filter elements 303 and 305 but not the filter element 
306. As a result, both the inverse amplifier and filter ele- 
ments 308 and 309 are used in the topology for the pre- 
5 distorter 300 of FIG. 11 but the inverse filter element 310 is 
not. 

Similar to the amplifier 301 of FIG. 9, the AM removal 
element 302 and the envelope detection element 304 of the 
pre -distorter 300 of FIG. 11 both receive an input signal 
10 x m(0 and split it into a PM (pulse modulation) signal x PA /t) 
and an envelope signal x e „ v (t). To compensate for the partial 
response f(x h (t) of the amplifier element 303, the corre- 
sponding inverse amplifier element 308 has a characteristic 
response D 1 (x(t)) to an input signal x(t). The output signal 
I 5 u(t) of the inverse amplifier element 308 in response to the 
envelope signal x env (t) is given by: 

u(t)=r 1 (x OT ,.(t)) (37) 

20 where x «,v(0= x (0- 

The inverse filter element 309 has a characteristic 
response h 1 _1 (t) to compensate for the response tq/t) of the 
filter element 305. The output signal v(t) of the inverse filter 
element 309 in response to the output signal u(t) from the 
25 inverse amplifier element 308 is given by: 

v(t)=hr 1 (t)*u(t)=hr 1 (t)*T 1 (x„ v (t)) (38) 

The mixer 307 then mixes the output signal u(t) from the 
inverse amplifier element 308 with the PM signal x PA /t) 
30 from the AM removal element 302. This produces the 
pre-distorted input signal x pred (t) of the pre-distorter 300. 

Similar to the manner in which the pre-distorters 100 and 
200 of FIGS. 5 and 8 were derived, the next step is to define 
the responses D 1 (x(t)) and h _1 (t) for the inverse amplifier 
35 and filter elements 308 and 309 of FIG. 11. In doing so, it 
is necessary to develop a model of the behavior of the 
corresponding amplifier and filter elements 303 and 305 of 
FIG. 9. The following discussion details how this is done. 

Like the filter element 102 of FIG. 4, the behavior of the 
40 filter element 305 can be modeled by an AR (i.e., auto 
regressive) filter if it is assumed to be minimum phase. Thus, 
the corresponding Z transform Hj(z) for the response h^t) 
can be modeled by the series expansion of Eq. (7), where 
Hj(z)=H(z) and the set of coefficients {a 0 , . . . ,a M } defines 
4:1 the response lq(t). This also means that the inverse response 
hj _1 (t) exists. The Z transform H 1 _1 (z) of this response 
hj _1 (t) of the inverse filter element 309 of the pre-distorter 
300 of FIG. 11 can be modeled by the series expansion given 
in Eq. (8). Thus, the set of coefficients {a 0 , . . . ,a M } also 
50 defines the response h 1 _1 (t). 

Since the response fq -1 )!) in fact exists, Eq. (33) can be 
re-written in the following manner: 

x„v(t)=x„„(t)*hr 1 (t) (39) 

55 

In view of Eq. (8), Eq. (39) can be re-written to model the 
behavior of the filter element 305 in an inverted manner 
according to: 

60 " ( 4 °) 

X env (t) S5 y, a m x con(t ~ m T) 
m = 0 

where 1/T represents the sampling rate. 

65 Eq. (35) can be re-written to model the behavior of the 
amplifier element 303 due to the partial response f(x fc (t)) 
using a series expansion as follows: 
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-continued 


,vW = fticontt)) ~ Yj bnXconVf 


(41) 


where the set of coefhcients {b 0 defines the partial 

response f(x fc (t)) of the amplifier element 303. Assuming 
that the inverse response f -1 (f(x co „(t))/g) to the response 
f( x con(t)) exists and is bounded by an input signal f(x co „(t))/ 
g, the behavior of the amplifier element 104 can also be 
modeled in an inverted manner with the following series 
expansion: 


xconlf) = r'OwW) 



(42) 


where the set of coefficients {d 0 , . . . ,d P } defines the 
response f' 1 (x(t)) of the inverse amplifier element 308 of the 
pre-distorter 300 of FIG. 11. The desired gain g of the 
amplifier element 303 appears in Eq. (42) so that the 
coefficients {d 0 , . . . ,d P } of the response T 1 (x(t)) are 
properly defined not to cause amplification of the pre- 
distorted output signal x pre /l) by the gain g. 

Using the assumption that the filter element 305 is mini- 
mum phase, the model of the amplifier 301 represented by 
Eqs. (40) to (42) is linear in its parameters. Thus, it can be 
easily and efficiently estimated using an LMS (i.e., least 
mean squares) estimation. 

In performing this LMS estimation, input matrices X conl , 
X con2 , and Y envl must first be defined as follows: 


XconiO X cim (t 7) ••• X con {I tllT) 
X C on (t + T) Xcon (0 Xcon (r-tm + llD) 

x C im(t + 27") x con (t + T) ••• x co „(t-(m+2)T ) 


(43) 


1 X con (t) ■■■ XcorVf 
1 X con (t + T) ■■■ X con (t + T) N 
1 x con (t + 2T) ■■■ Xcon (t + 2T) N 




1 tmv (r) 
t yem (r + T) 
1 y m ,(t + 2T) 


yenvi'f 
ytnv(l + T) N 
y m v(t + 2T) N 


(44) 


(45) 


Correspondingly, output vectors X env , Y env2 ) and X con3 
must then be defined in the following manner: 


X m 


XmM 
Xem(t + T) 
x env (r+ 2T) 


(46) 


V„v(0 


y r nM+T) 

v t „ v (r + 2T) 


(47) 


5 


Xco,,(i) 


XcoM + T) 

Xeoni* + 2T) 


(48) 


As will be discussed later in more detail, the elements of 
the matrices and vectors X con l , X con 2 , Y e „ v>1 , X em „ Y env>2 , 
1° and X con 3 are obtained by applying an input signal x, n (t) to 
the amplifier 301 of FIG. 9. Measurements of the signals 
y e „v> and x e „ v , y em ,' 2 are then made from which the 
elements of the matrices and vectors X conl , X con2 , Y envl , 
X env , Y mv2 , and X coll3 are obtained. These measurements 
15 are made at the sampling rate 1/T mentioned earlier. 

Eqs. (40) to (42) can now be written in matrix form 
according to the following relationships: 


20 


25 


30 


Xpnv ~ X con i 


’ env, 1 — Acon,2 


Xcon.i — y env, 2 


(49) 


(50) 


(5D 


In order to obtain the coefficient vectors [-a 0 . . . -a M Y, 
[b 0 . . . by] 7 ) and [d 0 . . . d p ] r , the LMS solution for Eqs. "(49) 
to (51) can be easily computed according to: 


35 


40 


45 


-a o 

-o« 

bo 


— ( X? 0 n,lXcon,l ) Xe on3 Xenv 


~ 1 Xe 0 n.lX C0n2 ) 1 X^^Yenv.l 


~ (^env,2^env,2) * X^nv.lXconA 


(52) 


(53) 


(54) 


where the matrices X con l H , X con2 H , and Y env2 H are the 
respective Hermitian transposes of the matrices X con l , X con 
50 2 , and Y env 2 . 

In this manner, the sets of coefficients {b 0 , . . . ,b^} and 
{a 0 , . . . ,& M } are estimated for properly defining the partial 
response f(x fc (t)) and hj(t) of the amplifier and filter elements 
303 and 305, respectively, of the amplifier 301 of FIG. 9. 
55 More importantly, the sets of coefficients {d 0 , . . . ,d p } and 
{a 0 , . . . ,& M } have been estimated for defining the responses 
f _1 (x(t)) and h 1 _1 (t) for the inverse amplifier and filter 
elements 308 and 309, respectively, of the pre-distorter 300 
of FIG. 11. Since the set of coefficients {d 0 , . . . ,d P } are 
60 directly estimated, the set of coefficients {b 0 , . . . ,b w } may 
be optionally estimated using Eq. (53) to determine the 
range over which the response f(x fc (t)) is invertible to obtain 
the inverse response f , (x(t)). 

The behavior of the inverse amplifier and filter elements 
65 308 and 309 can now be properly modeled by two cascaded 
processes. The first process is performed by the inverse 
amplifier element 308 and can be expressed as follows: 
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u(t ) K [ Jt mv (0° XoivU)' 



dp 


(55) 


The second process is performed by the inverse filter ele- 
ment 309 and can be expressed in the following manner: 


v(t)zi[u(t-T) u(f - 21) 


-ao 


u(t-MT)]\ 


-a i 


-om 


(56) 


Here as well, the processes only involve simple additions 
and multiplications. Thus, the implementation of the pre- 
distorter 300 with the inverse filter and amplifier elements 
308 and 309 is simple and efficient. 

Amplifier System 311 Topology and Operation 

FIG. 3 shows an amplifier system 311 that comprises the 
amplifier 301 and the corresponding pre-distorter 300 that 
was just derived. The pre-distorter 300 pre-distorts the 
original input signal x in (t) to produce the pre-distorted input 
signal x Jt) This done using the processes described in 
Eqs. (55) and (56). The pre-distorted input signal x JX) is 
then amplified by the amplifier 301 to produce the undis- 
torted output signal y ou /t) The output signal y OH ,(t) is not 
distorted by the non-linear behavior of the amplifier element 
303 and the bandlimiting effect of the filter element 305 that 
were compensated for with the pre-distorter 300. 

Performance Evaluation of Pre-Distorter 300 

Similar to the method used to derive the pre-distorter 100, 
the method just described for deriving the pre-distorter 300 
was verified with simulations. These simulations used a 
characteristic response f(x) for the non-linear amplifier 
element 303 of the amplifier 301 of FIG. 9 given by: 

f(x)=5x-2x 2 +0.01x 3 (57) 

The characteristic responses h 1 (t) and h 2 (t) of the filter 
elements 305 and 306 that were used had Z transforms Hj(z) 
and H 2 (z) given by: 


1 

Hl(i> ~ 1 - .09428 j- 1 + .3333i- 2 


HM = 


1 

l-.lsr 1 


(58) 

(59) 


In view of the responses f(x(t)) and h 1 (t) given by Eqs. 
(57) and (58), the filter element 305 will cause very little 
distortion while the amplifier element 303 will cause sig- 
nificant distortion if the input signal x,-„(t) is narrow band. 
On the other hand, if the input signal x,-,,)!) is wide band, the 
filter element 305 will cause significant distortion while the 
amplifier element 303 will not. 

The input signal x i;l (t) used for the simulations comprised 
a normalized two-tone single-sideband signal. To show the 
effects of the different types of distortion just described, the 
frequency spacing Af of the two tones was varied. 
Specifically, the frequency spacing Af varied between 0.01 
(1/T) and 0.10(f/T), where 1/T is the sampling rate. 

For each frequency spacing Af, a corresponding simula- 
tion of the amplifier 300 was performed using the responses 
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f(x(t)), hj(t), and h 2 (t) just described. Measurements of the 
signals x co „, y,,„ v , and x env , y e „ v>2 were made for each 
simulation for obtaining the elements of the matrices and 
vectors C con l , X con2 , Y envl , X env , Y e „ v ,, and X co „ j3 . Then, 
5 for each simulation, pre-distorters 300 with orders P=0 (i.e., 
no pre-distortion), 3, and 11 were derived in the manner 
described earlier. 

To evaluate the performances of the pre-distorters 300 for 
each simulation, the signal to distortion ratio given in Eq. 
to (26) is used. FIG. 12 shows these performances. When no 
pre -distortion is used (i.e., P=0), the SDR decreases as the 
signal frequency spacing Af increases. This is to be expected 
considering the bandlimiting effect of the response h 2 (t) of 
the filter element 305. For a relative frequency spacing Af of 
15 less than 0.04, the pre-distorter 300 of order P=ll is superior 
to that of order P=3. When the spacing is more that 0.04, the 
different pre-distorters 300 are equivalent due to frequency 
folding. 

20 Alternative Estimation Techniques 

The set of coefficients {d 0 , . . . ,d p } may be indirectly 
estimated in a similar manner to that described earlier for the 
pre -distorter 100 of FIG. 5. In this case, Eqs. (30) to (32) 
would also be used, except that the matrix Z 3 and the vector 
X 2 would be replaced by the matrix Y env2 and the vector 
X co „ i3 , respectively. 

Furthermore, as with the pre -distorter 100 of FIG. 5, RLS 
or GAL estimations may be used to obtain the coefficient 
30 vectors [b 0 . . . b^], [-a 0 . . . _ aM f, and [d 0 . . . d P ] r In this 
way, the matrix inversions used in the LMS estimations for 
solving Eqs. (52) to (54). Furthermore, SVDs (i.e., singular 
value decomposition) of the matrices X conl H X conl , 
K„„, 2 H Kon, 2 ’ and Yen V , 2 H Ye „ v ,2 ™ Eqs. (52) to (54) can be 
, 5 used for the same reasons as those expressed earlier for the 
pre-distorter 100 . 

Alternative Topology for Pre-Distorter 300 

In the earlier example, compensation was desired for the 
40 amplifier and filter elements 303 and 305 but not the filter 
element 306. As a result, corresponding inverse amplifier 
and filer elements 308 and 309 were utilized in the pre- 
distorter 300. 

However, as mentioned earlier, it may be desired to 
45 compensate for only one of or any combination of the 
amplifier element 303, the first filter element 305, and the 
second filter element 306. For example, if the bandlimiting 
effect of the filter element 306 also causes significant dis- 
tortion in the output signal y dis (t) at the carrier frequency, 
50 then compensation for it may be desired with the inverse 
filter element 310 of FIG. 10 and the pre-distorter 300 would 
have the topology shown in FIG. 11. In this case, the 
characteristic response h 2 _1 (t) of the inverse filter element 
310 would be derived in a similar manner to that in which 
55 the response h 1 _1 (t) was derived. 

Conclusion 

The method disclosed herein for deriving a corresponding 
60 pre-distorter for a non-linear amplifier has a number of 
distinct advantages. These are discussed in detail next. 

First, manufacturers do not have to make perfectly linear 
amplifiers because their non-linearities can be compensated 
for with corresponding pre-distorters. This allows the manu- 
65 facturers to make cheaper and simpler designs of amplifiers. 

Second, the method is adaptive since the method is 
dependent on making measurements of an amplifier to 
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derive the corresponding pre-distorter. Thus, the aging of the 
elements of an amplifier is not a problem since the pre- 
distorter can be continually adjusted based on new measure- 
ments. In this way, an amplifier system comprising the 
amplifier and the pre-distorter can stay linear over its life 5 
span. Similarly, the method can be applied to existing 
amplifiers in order to improve their linearity without any 
changes to the amplifiers. 

Finally, pre-distorters derived with this method are much 
easier to implement than those using the standard approach 10 
described earlier. Specifically, a pre-distorter derived using 
this method does not have to implement the kernels of the 
Volterra operator for the corresponding nonlinear amplifier 
and the inverse kernel of the first kernel of the Volterra 
operator for the non-linear amplifier. I 5 

While the present invention has been described with 
reference to a few specific embodiments, the description is 
illustrative of the invention and is not to be construed as 
limiting the invention. Various modifications may occur to 
those skilled in the art without departing from the true spirit 20 
and scope of the invention as defined by the appended 
claims. 

What is claimed is: 

1. A method of deriving a corresponding pre -distorter to 

a non-linear amplifier, the non-linear amplifier having a 25 
topology comprising filter and non-linear amplifier 
elements, the method comprising the steps of: 

defining the pre-distorter’s topology by: 

determining a corresponding inverse topology to the 
non-linear amplifier’s topology, the inverse topology 
comprising corresponding inverse filter and ampli- 
fier elements to the amplifier’s filter and non-linear 
amplifier elements; 

selecting each of the filter and non-linear amplifier 
elements for which compensation with the pre- 
distorter is desired; and 

removing from the inverse topology any of the inverse 
filter and amplifier elements corresponding to one of 
the filter and non-linear amplifier elements that was 4Q 
not identified so as to form the pre-distorter’s 
topology, the pre-distorter’s topology comprising 
each of the inverse filter and amplifier elements 
corresponding to one of the identified filter and 
non-linear elements; 

45 

define the identified inverse filter and amplifier elements’ 
responses by: 

applying an input signal to the non-linear amplifier and 
making measurements of selected signals of the 
non-linear amplifier; and 50 

estimating their respective sets of coefficients based on 
the measurements. 

2. The method of claim 1 wherein the input signal is the 

type of input signal that will be used in actual operating 
conditions for the non-linear amplifier. 55 

3. The method of claim 1 wherein at least one amplifier 
element is selected for compensation. 

4. The method of claim 1 wherein at least one filter 
element is selected for compensation. 

5. A corresponding pre-distorter to a non-linear amplifier, 60 
the non-linear amplifier having a topology comprising filter 
and non-linear amplifier elements, the pre-distorter compris- 
ing: 

a topology defined by: 

determining a corresponding inverse topology to the 65 
non-linear amplifier’s topology, the inverse topology 
comprising corresponding inverse filter and ampli- 
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fier elements to the amplifier’s filter and non-linear 
amplifier elements; 

selecting each of the filter and nonlinear amplifier 
elements for which compensation with the pre- 
distorter is desired; and 

removing from the inverse topology any of the inverse 
filter and amplifier elements corresponding to one of 
the filter and non-linear amplifier elements that was 
not identified so as to form the pre-distorter’s 
topology, the pre-distorter’s topology comprising 
each of the inverse filter and amplifier elements 
corresponding to one of the identified filter and 
non-linear elements; 

responses for the identified inverse filter and amplifier 
elements’ defined by: 

applying an input signal to the non-linear amplifier and 
making measurements of selected signals of the 
non-linear amplifier; and 

estimating their respective sets of coefficients based on 
the measurements. 

6 . The pre -distorter of claim 5 wherein the input signal is 
the type of input signal that will be used in actual operating 
conditions for the non-linear amplifier. 

7. The pre-distorter of claim 5 wherein at least one 
amplifier element is selected for compensation. 

8 . The pre-distorter of claim 5 wherein at least one filter 
element is selected for compensation. 

9. A method of compensating for distortion caused by a 
non-linear amplifier exhibiting filter and non-linear ampli- 
fier characteristics, the method comprising: 

deriving a pre-distorter, the pre-distorter exhibiting 
inverse filter and inverse non-linear amplifier 
characteristics, that acts on an input signal to generate 
a pre-distorted signal, wherein said pre-distorter com- 
pensates for distortion caused by the non-linear ampli- 
fier when said non-linear amplifier receives and ampli- 
fies said pre-distorted signal from said pre-distorter; 

defining the pre-distorter’s inverse filter and inverse non- 
linear amplifier characteristics by: 
developing a model of the non-linear amplifier’s filter 
and non-linear amplifier characteristics; 
developing a model of the pre-distorter’s correspond- 
ing inverse filter and inverse non-linear amplifier 
characteristics; 

applying a calibration signal directly to the non-linear 
amplifier and making measurements of selected sig- 
nals; and 

estimating coefficients for the pre-distorter’s inverse 
filter and inverse non-linear amplifier characteristics 
based on the measurements. 

10. The method of claim 9, wherein the non-linear ampli- 
fier’s first filter and non-linear amplifier characteristics cor- 
respond to filter and non-linear amplifier elements, and the 
pre-distorter’s filter and non-linear amplifier characteristics 
correspond to inverse filter and inverse non-linear amplifier 
elements, said pre-distorter’s elements and their arrange- 
ment being defined by: 

determining a corresponding inverse arrangement of ele- 
ments to the non-linear amplifier’s arrangements of 
elements; 

selecting each of the filter and non-linear amplifier ele- 
ments for which compensation with the pre-distorter is 
desired and; 

removing from the inverse arrangement any of the inverse 
filter and inverse non-linear amplifier elements corre- 
sponding to one of the filter and non-linear amplifier 
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elements that was not identified so as to form the 
pre-distorter’s arrangement of elements, the pre- 
distorter’s arrangement of elements comprising each of 
the inverse filter and inverse non-linear amplifier ele- 
ments corresponding to one of the selected filter and 5 
non-linear elements. 

11. The method of claim 9, wherein the calibration signal 
is the type of signal that will be used in actual operating 
conditions for the non-linear amplifier. 

12. The method of claim 10, wherein at least one amplifier to 
element is selected for compensation. 

13. The method of claim 10, wherein at least one filter 
element is selected for compensation. 

14. Apre-distorter to a non-linear amplifier, the non-linear 
amplifier having an arrangement of filter and non-linear 15 
amplifier elements, the pre-distorter comprising: 

an arrangement of inverse filter and inverse non-linear 
amplifier elements that acts on a signal to generate a 
pre-distorted signal, wherein said pre-distorter com- 
pensates for distortion caused by the non-linear ampli- 20 
fier when said non-linear amplifier amplifies said pre- 
distorted signal; 

the inverse filter and inverse non-linear amplifier elements 
having responses defined by: 

applying a calibration signal directly to the non-linear 25 
amplifier and making measurements of selected sig- 
nals; and 

estimating coefficients for the pre-distorter’s inverse 
filter and inverse non-linear amplifier elements based 
on the measurements. 30 

15. The pre-distorter of claim 14, wherein the pre- 
distorter comprises: 

inverse topology to the non-linear amplifier’s topology, 
the inverse topology comprising inverse filter and 
inverse non-linear amplifier elements corresponding to 
selected of the non-linear amplifier’s filter and non- 
linear amplifier elements; 

said inverse topology comprising an arrangement of said 
inverse filter and inverse non-linear amplifier elements 4(j 
opposite to the arrangement of filter and non-linear 
amplifier elements in the non-linear amplifier, such that 
the inverse filter and inverse non-linear amplifier ele- 
ments operate on a signal in the opposite order as the 
filter and non-linear amplifier elements to which they 4? 
correspond in the non-linear amplifier. 

16. The pre-distorter of claim 15, wherein the calibration 
signal is the type of signal that will be used in actual 
operating conditions for the non-linear amplifier. 

17. The pre-distorter of claim 15, wherein at least one 5Q 
amplifier element is selected for compensation. 

18. The pre-distorter of claim 15, wherein at least one 
filter element is selected for compensation. 

19. A pre-distorter for a non-linear amplifier, the non- 
linear amplifier comprising filter elements and non-linear 
amplifier elements, the pre-distorter comprising inverse fil- 
ter elements and inverse non-linear amplifier elements, 
wherein: 

the inverse filter elements are based on selected filter 
elements; 60 

the inverse non-linear amplifier elements are based on 
selected non-linear amplifier elements; 

each of the inverse filter elements correspond directly to 
one of the selected filter elements; 

each of the inverse non-linear amplifier elements corre- 65 
spond directly to one of the selected non-linear ampli- 
fier elements; 
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the inverse filter elements and the inverse non-linear 
amplifier elements are arranged within the pre-distorter 
such that they operate on a signal in opposite succes- 
sion as the selected filter elements and the selected 
non-linear amplifier elements to which they correspond 
within the non-linear amplifier; 
each inverse filter element exhibiting a response defined 
as the inverse of the response of the filter element to 
which it corresponds; and 

each inverse non-linear amplifier element exhibiting a 
response defined as the inverse of the response of the 
non-linear amplifier element to which it corresponds. 

20. The pre-distorter of claim 19, wherein: 

the selected filter elements and the selected non-linear 
amplifier elements are chosen based on their contribu- 
tion to distortion of a signal. 

21. The pre-distorter of claim 19, wherein: 

the inverse non-linear amplifier elements’ responses and 
the inverse filter elements’ responses are defined using 
measurements of selected signals after a calibration 
signal is received directly by the non-linear amplifier. 

22. The pre-distorter of claim 21, wherein: 
said definition is based, at least in part, on: 

models of the selected non-linear amplifier elements 
and selected filter elements based, at least in part, on 
a series representation of said elements; 
corresponding inverse models for the inverse non- 
linear amplifier elements and inverse filter elements 
based, at least in part, on a series representation of 
said elements; 

estimated coefficients for the inverse non-linear ampli- 
fier elements and inverse filter elements derived from 
said measurements. 

23. The pre-distorter of claim 22, wherein 

before said measurements are taken, choosing an order for 
the inverse non-linear amplifier elements’ models, the 
order determining how many coefficients are needed to 
describe the inverse non-linear amplifier elements. 

24. The pre-distorter of claim 21, wherein: 
the non-linear amplifier comprises: 

a filter element that receives a non-linear amplifier 
input signal and generates a first filtered signal; 
a non-linear amplifier element that receives the first 
filtered signal and generates a first output signal; 
the pre-distorter comprises: 

an inverse filter element that receives an input signal 
and generates an inverse filtered signal; 
an inverse non-linear amplifier element that receives an 
inverse filtered signal and generates the non-linear 
amplifier input signal. 

25. The pre-distorter of claim 24, wherein the measure- 
ments made comprise: 

samples of the calibration signal applied directly to the 
non-linear amplifier; and 
samples of the output signal. 

26. The pre-distorter of claim 21, wherein: 

the non-linear amplifier comprises two or more sub- 
amplifiers each sub-amplifier comprising: 
a filter element that receives a sub-amplifier input 
signal and generates a filtered signal; 
a non-linear amplifier element that receives the filtered 
signal and generates an output signal; 
the sub-amplifiers are cascaded to form the non-linear 
amplifier, wherein for each sub-amplifier, the output 
signal is the sub-amplifier input signal to the next 
sub-amplifier; 
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the pre-distorter comprises two or more corresponding 
sub-pre-distorters wherein each sub-pre-distorter com- 
prises: 

an inverse filter element that receives a sub-pre- 
distorter input signal and generates an inverse fil- 5 
tered signal; 

an inverse non-linear amplifier element that receives an 
inverse filtered signal and generates a sub-pre- 
distorter output signal; 

the sub-pre-distorters are then cascaded in opposite order 10 
to that of the sub-amplifiers to form the pre-distorter, 
for each sub-pre-distorter except final sub-pre-distorter, 
the sub-pre-distorter output signal is the sub-pre- 
distorter input signal to next sub-pre-distorter, the final 
output signal from the final sub-pre-distorter is the I s 
sub-amplifier input signal to be provided to the first 
sub-amplifier. 

27. The pre-distorter of claim 26, wherein the measure- 
ments made comprise: 

samples of the calibration signal applied directly to the 20 
non-linear amplifier; and 

samples of the output signal of each sub-amplifier. 

28. The pre-distorter of claim 21, wherein: 

the non-linear amplifier is an envelop elimination and 2 5 
restoration non-linear amplifier comprising: 
an AM removal element receiving a non-linear ampli- 
fier input signal and generating a pulse-modulated 
signal; 

an envelope detection element receiving a non-linear 30 
amplifier input signal and generating an envelope 
signal; 

a first filter element receiving an envelope signal and 
generating a control signal; 

a non-linear amplifier element receiving the pulse- 35 
modulated signal and the control signal and gener- 
ating an amplified signal; 

a second filter element receiving the amplified signal 
and generating a non-linear amplifier output signal; 

the pre-distorter comprising: 40 

an AM removal element receiving a pre-distorter input 
signal and generating a pulse-modulated signal; 
an envelope detection element receiving a pre-distorter 
input signal and generating an envelope signal; 
an inverse non-linear amplifier element receiving the 45 
envelope signal and generating a corrected envelope 
signal; 

a first inverse filter element receiving the corrected 
envelope signal and generating a filtered corrected 
envelope signal; 50 

a mixer receiving the pulse-modulated signal and the 
filtered corrected envelope signal and generating the 
non-linear amplifier input signals. 

29. The pre-distorter of claim 28, the pre-distorter addi- 
tionally comprising: 55 

a second inverse filter element receiving an input signal 
and generating the pre-distorter input signal. 

30. The pre-distorter of claim 28, the measurements 
comprising: 

samples of the calibration signal applied directly to the 60 
non-linear amplifier; 

samples of the control signal; 

samples of the envelope signal; 

samples of the amplified signal. 65 

31. A method for deriving a pre-distorter for a non-linear 
amplifier, the non-linear amplifier comprising filter and 
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non-linear amplifier elements, the pre-distorter comprising 
inverse filter and inverse non-linear amplifier elements, by: 

determining a corresponding inverse arrangement of ele- 
ments to the non-linear amplifier’s arrangement of 
elements, the inverse arrangement comprising corre- 
sponding inverse filter and inverse non-linear amplifier 
elements to the amplifier’s filter and non-linear ampli- 
fier elements, said inverse filter and inverse non-linear 
elements arranged in reverse order to their correspond- 
ing filter and non-linear amplifier elements, such that 
the last element to operate on a signal in the pre- 
distorter corresponds to the first element to operate on 
a signal in the non-linear amplifier; 

selecting each of the filter and non-linear amplifier ele- 
ments for which compensation with the pre-distorter is 
desired; 

removing from the inverse arrangement of elements any 
of the inverse filter and inverse non-linear amplifier 
elements corresponding to one of the filter and non- 
linear amplifier elements that was not selected, such 
that the pre-distorter comprises inverse filter and 
inverse non-linear amplifier elements corresponding 
only to selected filter and non-linear amplifier 
elements, 

defining the inverse filter and inverse non-linear amplifier 
elements’ responses by: 

applying a calibration signal directly to the non-linear 
amplifier and making measurements of selected sig- 
nals; and 

estimating their respective sets of coefficients based on 
the measurements. 

32. The method of claim 31, wherein: 

the selected filter and non-linear amplifier elements are 
selected based on their contribution to the distortion of 
an input signal. 

33. The method of claim 31, wherein: 

defining the inverse filter and inverse non-linear amplifier 
elements’ responses also comprises: 
developing a model for the selected non-linear ampli- 
fier’s filter and non-linear amplifier elements based, 
at least in part, on a series representation; 
developing a model for the inverse filter and inverse 
non-linear amplifier elements based, at least in part, 
on a series representation; and 
choosing an order for the inverse non-linear amplifier 
elements, the order determining the number of coef- 
ficients necessary to describe the non-linear ampli- 
fier elements’ responses. 

34. The method of claim 31, wherein: 

the calibration signal is based, at least in part, on the 
signal that will be used in actual operating conditions 
for the non-linear amplifier. 

35. The method of claim 31, wherein: 

estimating the coefficients is based, at least in part, on a 
least mean squares estimation. 

36. The method of claim 31, wherein: 

estimating the coefficients is based, at least in part, on a 
recursive least squares estimation. 

37. The method of claim 31, wherein: 

estimating the coefficients is based, at least in part, on a 
gradient adaptive lattice estimation. 

38. The method of claim 31, wherein: 

estimating the coefficients involves utilizing singular 
value decompositions. 





